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The boundary of two mixed Bose-Einstein condensates 
interacting repulsively was considered in the case of spatial 
separation at zero temperature. Analytical expressions for 
density distribution of condensates were obtained by solving 
two coupled nonlinear Gross-Pitaevskii equations in cases cor- 
responding weak and strong separation. These expressions 
allow to consider excitation spectrum of a particle confined in 
the vicinity of the boundary as well as surface waves associ- 
ated with surface tension. 

PACS Number(s): 03.75.Fi 



I. INTRODUCTION 

The experimental realization of Bose Einstein conden- 
sation in trapped dilute gases has allowed to inves- 
tigate variety of properties of quantum fluids both theo- 
retically and experimentally. In the recent years, it has 
become possible to produce and explore mixtures of Bose- 
Einstein condensates corresponding to different internal 
states 

Theoretical treatment of mixtures ]7|-|lC| assures that 
depending on the relative strength of interactions inside 
each condensate and between them it is possible to ob- 
serve spatial separation. Experimental realization of such 
systems [^-^] has given an opportunity to study both 
equilibrium properties and dynamics of separation. Al- 
though, structure of the boundary has been qualitatively 
analyzed in an analytical distribution of densities 
has not been derived yet. Authors of || discussed asymp- 
totic behavior of densities far from the boundary, and on 
the basis of these estimations analyzed surface tension. 

In the present paper, we explore the boundary between 
two repulsively interacting condensates at zero temper- 
ature in two limit cases corresponding weak and strong 
separation of the condensates. We show that in the case 
of weak separation it is possible to derive equations for 
densities, and using iteration method solve them analyt- 
ically. We show that asymptotic behavior of our solu- 
tion coincides with those predicted in [||. For the sake 
of completness, we also provide the solutions of these 
equations in the case of strong separation. Importance 
of the solution lies in the possibility to explore quanti- 
tatively different types of excitations on the boundary. 
The structure of the boundary in both cases allows to 
consider one-particle excitations as well as surface waves 
associated with the boundary. Obtained expressions for 
dispersion relation of surface waves can be used to ex- 
plore that phenomenon experimentally. 

The Hamiltonian describing the mixture of two weakly 



interacting Bose-Einstein condensates can be written in 
the form 
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Here, Uj = Airh ai/rrii > characterizes the inter- 
action inside each condensate, ui 2 — 1 r Kh 1 a\ 2 (m\ + 
TO 2 )/(mim 2 ) > - the intercondensate interaction, 
- mass of a particle of each condensate, a%, an - corre- 
sponding scattering lengths, Vi (r) - external trapping po- 
tentials. Theoretical treatment of the mixtures |7|,[| has 
shown that separation takes place when u\ 2 j ^Ju\u 2 > 1. 

Starting with the Hamiltonian ([!]) we get Gross- 
Pitaevskii equations for condensate wave functions: 
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As we are interested in studying of stationary solu- 
tions of these equations, then assuming as usual tyj cx 
exp(— ifj,jt), where fij - chemical potentials of the con- 
densates, we obtain two coupled nonlinear equations for 
densities of gases n»(r) = |^j(r)| 2 : 
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These equations are essentially nonlinear, so to find so- 
lutions we need to make some simplifications. We assume 
that Vi(r) = Va(r) and consider the case when the size 
of a boundary between condensates much less than char- 
acteristic length of the trap. In the case of a parabolic 
trap potential and Thomas-Fermi regime, it means that 
d <C Rt f , where d is the size of the boundary, and Rtf is 
the Thomas- Fermi radius of the atomic cloud. Physically, 
it helps to avoid the effect of the potential on the form of 
the boundary. To simplify calculations further, we also 
suppose that separation takes place in one dimension, z. 
This allows us to write the system of equations (^) in the 
following form: 
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Although, there is no explicit trapping potential in 
these equations, we have to impose external conditions on 
the solutions, that implicitly take it into account. Let us 
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assume that separation takes place along the z direction, 
and the condensate with the label "1" is to the right, and 
the one with the label "2" is to the left of the boundary. 
Then, asymptotically we require: 
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where mo, mo ~ equilibrium densities of condensates 
far from the boundary. Substituting these conditions in 
equations ^ for densities, and using the obvious condi- 
tion of equilibrium we obtain 



(j,! = uinio, P2 = u 2 n 20 , 

Pi = uinf /2 = P 2 = u 2 nl /2 



(6) 



Here, we used well-known expression for the pressure of 
a homogeneous weakly-interacting Bose gas. That con- 
dition connects equilibrium densities of condensates far 
from the boundary. 

To reduce the number of parameters in equations (0) 
we notice that it is possible to exclude difference of 
masses by the change: 



uimi/m 2 , u 2 = u 2 m 2 /rrii 



niy/m 2 /m 1 , n 2 = n 2 a/toi/to 2 , 
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We get the same equations (^) and conditions (||), (||) 
but for quantities with stars and with the same mass 
to*. To simplify notations, in the following discussion we 
omit stars. The generalization for different masses can 
be easily done by following the above rules. 

We can solve equations for densities (^) with condi- 
tions (|B|) , (^|) analytically in two limit cases for weak sep- 
aration when A = u\ 2 j ^Ju\u 2 — 1 <C 1, and strong sep- 
aration when A ^> 1. Notice, that both parameters u% 2 
and A are not affected by the above procedure of mass 
difference excluding. 



Densities of condensates "1" and "2" are easily ob- 
tained if functions p and g are known: 



»i = (w 2 M) 1/4 p[l+5]/2, 
n 2 = {u 1 /u 2 ) l ' A p[l-g\/2 
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It is straightforward to derive equations for p and g: 
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Here, /' = £of^, and we also introduced a ■ 

and = m{uiu 2 ) 1 / i {^/ul+ y^)p /h 2 . 

To find asymptotic solutions of equations (|ll|) in the 
case of A C 1 we use iteration method. Let us suppose 
that terms with derivatives of p are much smaller than 
the others. We will justify this assumption in the end of 
calculations. 

Neglecting terms with derivatives of p, from equa- 
tions (||) we get 
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The equation for g can be solved by substitution g' = 
/, so taking into account conditions (|9j) we get that g is 
the solution of the equation 
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At this point it is clear that assumptions we made were 
correct. Namely, we obtain that ^fp' cx A 3 / 2 , y/p" oc A 2 , 
g' oc •y/A, so neglected terms by A times less than the 
others. 

Now we can write down solutions for both p and g in 
the parametric form 



II. WEAK SEPARATION 

Let us consider weak separation when condition A = 
u 12 j ^Ju\u 2 — 1 <C 1 is satisfied. We expect that in the 
simplest case when u\ = u 2 the whole density of a gas 
n(z) = ni(z) + 712(2) is approximately constant. That is 
why, to find the solution it is natural instead of n± , n 2 to 
introduce other quantities and solve equations (Q) using 
small parameter A. Consider the functions: 



p = (wi/u 2 ) 1/4 ni + (u 2 /«i) 1/4 n 2 , 
g = [(ui/u 2 ) 1/4 ni - (u 2 /ui) 1/4 ri 2 ] jp 



(8) 



Conditions (||), (^) give us simple asymptotic behavior 
of these functions: 
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Here, the second equation is the solution of (|l3|). In these 
expressions an arbitrary constant z defines the position 
of the boundary, and in the case of finite geometry with 
given average number of particles in condensates can be 
obtained from the condition of equal pressures (^]) . There 
is also physically obvious symmetry z — zq — > zq — z, 

.9 -> -9, ot -> -a- 

Finally, the densities of condensates have the following 
form: 
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With the parametric equation for <?, these densities are 
the main result of the paper. On fig With the use of the 
desribed method, it is possible to derive expressions for 
densities in the next orders of the small parameter A in 
the form of asymptotic series. As follows from the above 
estimations of neglected terms in equations ([n]) the next 
order is proportional to A 2 . The typical dependence of 
densities on the distance from the boundary is shown in 
Fig. g and Fig. §. 

It is interesting to notice that the whole density n = 
?ii + n 2 at some values of parameters has a well on the 
boundary as shown in Fig. [|. In the case of weak separa- 
tion we see that it is broad and shallow. This result is a 
consequence of the interparticle interactions. The inter- 
particle potential for each condensate acts in some sense 
as a wall, so we expect that probability of a particle to 
be close to the boundary decreases, which means lower 
density. The existence of such a well allows to consider 
the possibility of confining of a particle of another sort in 
the vicinity of the boundary, which is discussed below. 

Although, we can not further simplify solutions (14), 
it is interesting to derive asymptotic behavior of these 
functions in particular limits: 
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Here, = %/ y/2rrii^i - correlation lengths of conden- 
sates, defined by chemical potential fXi and /i 2 , and we 
used (0) to include mass difference. The asymptotic be- 
havior far from the boundary of each condensate is de- 
fined by its correlation length. As we see, the size of the 
boundary can be approximated as d ~ (£i +^2)/VA and 
for A -C 1 appears to be much larger than correlation 
lengths. 

Solutions have the simplest form when a — 0: 
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Here, we choose z = 0. 



III. STRONG SEPARATION 

To analyze the case of strong separation A = 
U12/ y/u\U2 3> 1 (we use the same notation but for an- 
other quantity) we start from density equations @). In 
this case we expect that density on the boundary will 
be approximately zero because interparticle interactions 
make it almost impossible for one condensate to pene- 
trate inside the other. To estimate the density of conden- 
sates on the boundary we can use the fact that second 
derivatives of wave functions should be approximately 



zero there. That makes system 
equations with solution: 
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This allows us in the zero approximation to use simple 
conditions for the densities ri\{z < 0) = n 2 (z > 0) = 0. 
Then equations M) have simple form: 
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Solutions are easily obtained: 



ni{z > 0) = ni tanh 2 [-^; 
n 2 {z < 0) = n 20 tanh 2 [^-' 



(20) 



Here, nio and n 2 o are connected by the condition of equal 
pressures @ , and we choose the position of the boundary 
at z = 0. 

The size of the boundary in this case is approximately 
d w 2\/2(£i + £2)- The dependence of the densities on 
the distance from the boundary is shown in Fig. ^[ As 
in the previous section, we see that there is again a well 
in the whole density but in the case of strong separation 
it becomes narrower and deeper in comparison with the 
one for weak separation. 



IV. ONE-PARTICLE EXCITATIONS ON THE 
BOUNDARY 

The existence of a well in the whole density allows to 
consider confining of a particle of another sort in the 
vicinity of the boundary. As a general property of a 
quantum-mechanical motion in a one-dimensional well 
there always exists such a confined state. As an exam- 
ple, we consider the simplest case when a = and a 
particle of another sort interacts with both condensates 
repulsively with the same constant A. The Schrodinger 
equation for the wave function of a particle with the mass 
M has the form: 
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We can solve equation ( |2l| ) for weak and strong sepa- 
ration cases simultaneously. It has universal form 
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where e = E — Xuq, and Uq = 3AAno/4, f3 = vA/£o for 
weak; and Uq = Arto, ft = l/(\/2£o) for strong separation. 
The spectrum of energy e is well-known: 



3 



-u— 

^ 4AI 



-n 



8.1/ 



>1- 


|_ g MA 




mu 







"weak" , 
2 (23) 
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where p = utiq - chemical potential, j = 0, 1, . . . and the 
condition that an expression in [. . .] is positive defines 
the upper limit for j. There is always at least one state 
with j = 0. 



V. SURFACE WAVES ON THE BOUNDARY 

There is also another type of excitations associated 
with the boundary. As we see condensate density dis- 
tributions give rise to nonzero surface tension which was 
previously analyzed in Q , and we use the expression for 
surface tension derived there: 



-t-oo 



dz 
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Substituting expressions (^5|) for densities in the case 
of weak separation and taking into account only first 
nonzero order in A we obtain for surface tension 



6 



, 1-yT^Z? 



(25) 



where a = (mi^/u{ — rti2^Ju2) / ( m iy/ui + m i\f^ : i) an(1 
we used ^ to include mass difference; £i, £2 - correlation 
lengths of condensates; P - the pressure given by (|^). 
When a = 0, a w = Pv^A(& + &)/4 

In the case of strong separation we use expressions ( poj ) 
to get surface tension 



PV2 
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Let us notice that expressions (|25|) and (|26|) differ 
from the ones obtained in B. Although, in qualitative 
sense our expressions coincide with those of S , using our 
method of solution we can get a general expression appli- 
cable for variety of parameters, and for example retrieve 
the correct numerical factor for the case considered in || . 
As follows from the general expression for surface ten- 
sion (24), we need to know the behavior of densities not 
only far but also in the vicinity of the boundary. That is 
why, estimate character of expressions for densities in ||] 
could only give qualitative answer for surface tension. 

For velocities smaller than the speed of sound we can 
consider gas as incompressible, so it is possible to write 
down the usual hydrodynamic equations and find disper- 
sion relation of surface waves on the boundary associated 
with surface tension. Suppose that we have condensates 
in a "box" and the position of the boundary is defined 
by the distances L\ and L2 from the "box" walls, where 



labels "1" , "2" correspond to the side with the same con- 
densate, and the box sizes for other directions are L x , 
L y . Taking into account the fact that velocities are zero 
on the walls we get the dispersion relation 

ak 3 
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Pi coth(fcLi) + P2 coth(/cL2) ' 
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where a is the surface tension; pio = minio, P20 = rn2U2o 
- mass densities; wave vector along the boundary k = 
^f(irn x /L x ) 2 + (Trriy/Ly) 2 , n x ,n y = 0,1..., not n x = 
n y = 0; Li, L2 - sizes of condensates in the direction 
perpendicular to the boundary. We assume that Li, L2 
are much larger than the size of the boundary d. 

For dispersion relation we consider two limit cases cor- 
responding to wavelengths A < L, (fcij 3> 1) and A ^> L, 
(kLi <C 1). Let us notice that the second case is possible 
only if Li < L XiV . 

In the first case we obtain 
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In the long wavelength limit we get 
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As we see, for weak separation uj oc A 1 / 4 , and surface 
waves are relatively "soft" in that case, which enables to 
consider them as a dissipative channel for other conden- 
sate excitations. 



VI. CONCLUSION 

We performed the analysis of the boundary of two 
Bose-Einstein condensates interacting repulsively in the 
limit cases corresponding weak and strong separation at 
zero temperature. 

For weak separation we obtained solutions (|l^) of two 
coupled nonlinear Gross-Pitaevskii equations using small 
parameter A. The solutions show that the penetration 
depth of the condensate "i" inside the other is estimated 
as fi/VA, so that the size of the boundary d ~ (£1 + 
£2)/VA is much larger than correlation lengths, which 
was obtained experimentally [pi. There is also a well in 
the full density profile, which is the consequence of the 
wave function behavior near the boundary. In general, 
the proposed method of obtaining density distributions 
for the weak separation case if desired can be expanded 
to obtain expressions for the next orders of A. 

We also considered the case of strong condensate sep- 
aration but restricted ourselves to the zero order ap- 
proximation. In that case the size of the boundary is 
d ~ (£1 + £2), and the whole density of gases goes ap- 
proximately to zero on the boundary. 
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The existence of the well in the density profile at some 
parameters allowed us to consider one-particle excita- 
tions on the boundary. Using the expressions for den- 
sities we found excitation spectrum of a particle in the 
simplest case when constants of interaction of the par- 
ticle with condensates are the same, and distribution of 
densities on the boundary corresponds to a = 0. The 
generalization to other cases can be easily done with the 
use of expressions (|l5]) , (|2C]) for density distributions. Let 
us notice that the existence of a potential well depends on 
the interaction constants of a particle with condensates 
as well as the relation between interaction constants of 
the condensates. Observation of such confined states can 
be possible only if temperature are smaller than the po- 
tential well. 

It was shown that there also exist collective excita- 
tions associated with surface tension. The expressions for 
surface tension were obtained for both weak and strong 
separation. The dispersion relation for surface waves was 
analysed in the case where condensates fill finite volumes. 
The dispersion relation has different forms in two cases 
corresponding to short- and long- wavelength limit. In 
the case of weak separation "soft" surface modes can be 
a dissipative channel of other condensate excitations. 

Author gratefully acknowledges valuable discussions 
with Yu.M. Kagan and other members of the Theoret- 
ical Division of RRC "Kurchatov Institute" , Russia. 
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FIG. 1. Distribution of densities in the case of weak sep- 
aration. A = 0.3, a — 0.1. Densities are in the scale of 
nio. Distance is in the scale of fo/-\/A(l — a 2 ). There is a 
well in the whole density profile. 
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FIG. 2. Distribution of densities in the case of weak sep- 
aration. A = 0.3, a = 0.3. Densities are in the scale of nio. 
Distance is in the scale of £ /^A(l - a 2 ) 
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FIG. 3. Distribution of densities in the case of strong 
separation when ui/«2 = 4. Densities are in the scale of nio- 
Distance is in the scale of \/2£i. 
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